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ABSTRACT 

CyclotrimethlenetrinitraimLne  (RDX  or  Cyclonite)  exists  in  two  poly- 
morphic species,  a  and  3.  The  a-form  is  stable  under  normal  conditions, 
whereas  3-RDX  is  formed  by  recrystallization  from  high  boiling  solvents. 
It  immediately  transfers  to  the  a-form  in  the  presence  of  excess  a-RDX, 
and  it  is  therefore  assumed  that  RDX  as  it  exists  in  its  normal  state 
consists  solely  of  the  a  polymorph. 

The  thermal  decomposition  of  RDX  has  been  comprehensively  studied 
at  Picatinny  Arsenal  by  means  of  quantitative  measurements  of  the  de- 
composition products.  The  thermal  stability  of  explosive  compositions 
has  also  been  investigated  by  measuring  the  rate  of  gas  evolution  during 
decomposition. 

It  was  the  purpose  of  this  work  to  obtain  thermal  and  kinetic  data 
for  RDX  using  Differential  Thermal  Analysis  (DTA)  experimental  methods. 
These  results  and  techniques  were  compared  with  those  obtained  using 
the  other  methods  and  with  theoretical  values.  The  favorable  results 
lend  credibility  to  the  DTA  method  and  suggest  that  further  development 
in  the  technique  will  prove  useful  in  providing  a  fast  and  reasonably 
accurate  method  of  characterizing  the  thermal  and  kinetic  properties  of 
new  explosives  and  propellants. 
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I.  INTRODUCTION 

Explosives  undergo  thermal  decomposition  at  temperatures  far  below 
those  at  which  spontaneous  explosion  may  occur.  The  decomposition  re- 
action is  important  in  determining  the  stability  of  the  explosive.  A 
considerable  amount  of  heat  energy  is  liberated  in  the  decomposition 
process.  This  highly  exothermic  reaction  accelerates  to  a  high  veloc- 
ity culminating  in  deflagration  or  detonation  when  the  rate  of  heat 
production  is  greater  than  the  rate  of  heat  loss  to  the  surroundings 
by  conduction  and  convection.  Explosives  in  their  normal  state  are 
metastable,  the  free  energy  of  the  explosive  being  considerably  higher 
than  its  decomposition  products.  For  an  explosive  to  remain  stable  at 
ordinary  temperatures,  its  decomposition  reaction  must  depend  on  rel- 
atively high  activation  energies. 

Kinetic  studies  are  therefore  important  in  the  study  of  explosive 
reaction  mechanisms.  The  Differential  Thermal  Analysis  technique  is 
one  of  several  that  is  used  to  study  reaction  kinetics.  When  a  reac- 
tion occurs  in  DTA,  the  change  in  heat  content  and  in  the  thermal  prop- 
erties of  the  sample  is  indicated  by  a  deflection  in  the  thermogram. 
For  reactions  that  possess  activation  energies  and  thereby  proceed  at 
a  rate  varying  with  temperature,  the  position  of  the  peak  varies  with 
the  heating  rate.  The  activation  energy  and  Arrhenius  frequency  factor 
can  be  determined  from  this  experimental  information. 

Cyclotrimethylenetrinitramine  (RDX  or  Cyclonite)  is  an  important  ex- 
plosive used  in  many  military  composite  explosives.  RDX  melts  at  203°C, 
and  the  decomposition  of  liquid  cyclonite  follows  the  unimolecular  rate 


equation  over  the  range  of  temperatures  213-299°C   [Robertson  19*19]  • 
The  Arrhenius  frequency  factor  is  large,  and  the  activation  energy  is 
of  the  order  of  50  kcal/mole.  The  high  value  of  the  Arrhenius  frequency 
factor  is  a  possible  indication  of  self  heating. 

The  purpose  of  this  work  was  to  investigate  the  thermal  decomposi- 
tion of  RDX  using  DTA  experimental  techniques.  In  the  decomposition 
of  lj3j5j7-Tetranitro-l,3i5,7-Tetrazacyclooctane  (HMX),  three  separate 
peaks  were  observed  in  the  DTA  thermogram,  and  three  separate  activa- 
tion energies  obtained  from  the  DTA  data  [Hondee  1971].  It  was  of 
interest  to  determine  whether  these  phenomena  were  evident  in  the  de- 
composition of  RDX,  to  detect  any  evidence  of  self-heating,  and  to 
compare  the  kinetic  data  obtained  with  that  found  in  the  literature. 


II.  THERMAL  DECOMPOSITION  OF  EXPLOSIVES 

The  investigation  of  the  thermal  behavior  of  explosives  involves  both 

the  study  of  chemical  kinetics  and  the  theory  of  transport  properties. 

The  chemical  kinetics  of  explosive  reactions  are  not  fully  understood  due 

to  the  high  reaction  rates  and  temperatures  involved  and  the  complexity 

of  the  products  of  the  reaction.  The  decomposition  of  RDX  may  be  written: 

C  H  0  N  +  0  ■*-  3C0  +  3H  0  +  6N0 
3666      2        2       2  2 

This  reaction  may  be  accurate  for  detonation  but  is  unlikely  to  occur 
thermally  due  to  the  large  number  of  bonds  which  must  be  broken  in  the 
parent  molecule.  Attempts  to  study  these  reactions  require  that  mea- 
surements of  the  decomposition  rates  be  made  at  relatively  low 
temperatures . 

Thermal  decomposition  is  essentially  a  heat-balance  problem.  The 
heat  loss  to  the  surroundings  and  the  accumulation  of  heat  in  the  ex- 
plosive equals  the  chemical  energy  generated  by  the  decomposition  of 

the  explosive. 

<V2T  +  pc<j>  =  dq/dt 
This  can  be  written  simply  as  F  +  G  =  H.  The  chemical  energy  generated 
by  the  decomposition,  H,  increases  exponentially  with  temperature.  The 
heat  loss,  F,  consists  of  thermal  conduction,  which  is  proportional  to 
the  thermal  gradient,  and  radiation,  which  follows  the  Stefan-Bolt zmann 
T1*  law.  The  heat  loss  therefore  increases  at  a  slower  rate,  especially 
since  the  heats  of  activation  of  explosives  are  of  the  order  of  50 
kcal/mole.  The  accumulation  of  heat  in  the  explosive,  G,  becomes  more 
important  as  the  temperature  increases,  causing  acceleration  of  the 
decomposition. 
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Some  explosives  liberate  energy  in  their  interior  due  to  a  slow  chem- 
ical reaction  that  is  often  a  decomposition  reaction.  This  is  known  as 
self-heating  and  causes  the  temperature  of  the  material  to  increase,  re- 
sulting in  an  increase  in  the  rate  of  reaction.  Dynamic  equilibrium  is 
reached  when  the  heat  is  removed  as  fast  as  it  is  generated.  If  heat  is 
generated  faster  than  it  can  be  removed,  the  temperature  increases  more 
rapidly  and  the  reaction  accelerates.  The  maximum  temperature  of  the 
surroundings  for  which  a  steady  state  is  possible  is  called  the  critical 
temperature  for  the  explosive.  Above  this  temperature,  the  material  will 
self-heat  until  ignition  or  explosion  occurs  [Longwell  196I]. 

Decomposition  reactions  are  sometimes  autocatalytic  or  catalyzed  by 
small  amounts  of  impurities.  The  decomposition  products  or  impurities 
in  the  explosive  may  act  as  a  catalyst  in  that  the  reaction  rate  is  in- 
creased as  more  products  are  produced.  The  decomposition  rates  of  these 
materials  then  become  a  function  of  the  duration  or  extent  of  the 
decomposition. 
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III.   REACTION  KINETICS  IN  DTA 

The  material  of  this  section  is  a  summary  of  Kissinger's  work  on 
the  theory  of  differential  thermal  analysis  [Ref.  2].  DTA  techniques 
produce  a  deflection  or  peak  (AT)  in  the  thermogram  due  to  a  change  in 
heat  content  and  thermal  properties  of  the  active  sample  during  a  re- 
action. This  peak  temperature  is  affected  by  experimental  technique, 
and  is  often  higher  than  known  transition  or  decomposition  temperatures. 

The  temperature  of  maximum  deflection  is  defined  by  kinetic  con- 
stants and  the  heating  rate.  If  a  reaction  posses  an  activation  energy, 
it  will  proceed  at  a  rate  varying  with  temperature  and  the  position  of 
the  peak  will  vary  with  the  heating  rate.  The  temperature  of  maximum 
deflection  in  DTA  is  also  the  temperature  at  which  the  reaction  rate  is 
a  maximum.  The  temperature  distribution  of  the  sample  explosive  obeys 
the  general  heat  flow  equation; 

3t   pc     "  pc  dt  U^ 

where  the  rate  of  heat  generated  due  to  chemical  reaction  per  unit  vol- 
ume of  the  sample  is  equal  to  zero  in  the  inert  reference.  The  temper- 
ature distribution  in  the  reference  material  is  then  given  by; 

H  -  -r v2T  <2) 

9t   pc 
Assuming  the  sample  to  be  a  cylinder  of  radius  a  and  of  infinite 
length,  integration  of  equation  2  yields; 

Tr  =  TQ  +  ft  _^§i  (3) 

where  <f>  =  dT/dt  and  the  outside  temperature  is  given  by  T  =  T0  +  <j>t. 
The  rate  of  heat  generated  due  to  chemical  reaction  is  a  function  of 
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temperature  in  the  active  sample.  Equation  1  is  therefore  a  non-linear 
partial  differential  equation.  Assuming  that  the  temperature  of  the  out- 
side of  the  tube  rises  at  a  linear  rate,  the  solution  for  the  temperature 
at  the  center  of  the  sample  will  be  of  the  form; 

Ts  =  Th  +  4>t  -  f(dq/dt)  (4) 

where  f(dq/dt)  is  a  function  of  the  reaction  rate  and  includes  secondary 
effects  of  the  reaction  such  as  changes  in  volume,  density  and  thermal 
properties . 

The  differential  temperature  is  the  difference  between  temperatures 
at  the  centers  of  the  two  samples; 

fl  -  ffe9\      _  /4>pca2\ 

D    \dt /sample   \  k  /reference  (5) 

and  is  a  maximum  at  the.  temperature  of  maximum  deflection  where  d6/dt 
is  zero. 

de/dt  =  f  (dq/dt)  d2q/dt2  (6) 

When  the  derivative  of  the  rate  of  heat  absorption  is  zero  it  can  be 
seen  from  Equation  6  that  de/dt  is  also  zero.  Since  the  rate  of  heat 
absorption  is  proportional  to  the  rate  of  reaction,  Equation  6  shows 
that  the  peak  differential  deflection  occurs  when  the  reaction  rate  is 
a  maximum.  It  should  be  noted  that  the  heating  rate  must  be  constant 
for  this  to  hold  and  that  Equation  6  is  valid  for  a  sample  of  any  shape. 

The  rate  law  for  the  thermal  decomposition  of  explosives  can  be 
described  by  the  equation; 

Ox/3t)T  =  kp  (l-x)n  (7) 

where  fcj  =  A  exp(-Ea/RT)  and  temperature  is  a  function  of  time. 

^L=  f^\     +  (*A      dT  (8) 

dt   ^3t/T   \3T/t  dt 

13 


The  rate  of  change  of  x  with  respect  to  temperature,  with  the  time 
coordinates  fixed,  is  equal  to  zero  because  the  number  and  positions 
of  the  particles  of  the  reactant  are  also  fixed  with  time.  The  only 
effect  of  instantaneous  change  in  temperature  is  in  the  velocities  of 
thermal  motion  of  the  reactant  particles.  Equation  7  can  now  be 
written; 

dx/dt  -  A(l-x)nexp(-Ea/RT)  (9) 

which  holds  for  any  temperature  as  long  as  the  fraction  of  reactant 
decomposed  and  temperature  are  measured  at  the  same  instant.  It  is 
assumed  that  the  empirical  order  of  the  reaction  remains  constant 
through  the  greater  part  of  the  reaction. 

when  the  temperature  rises  during  the  reaction,  the  reaction  rate 
will  rise  to  a  maximum 'value,  then  return  to  zero  as  the  reactant  is 
exhausted.  This  maximum  rate  occurs  when  d/dt (dx/dt)  is  zero.  Pro- 
viding the  temperature  rises  at  a  constant  rate,  the  differentiation 
of  Equation  9  yields; 


d_/dx\   dx 
dt\dtj  "  dt 


M   -  An(l-x)n-l  e-Ea/RT 

RT 


(10) 

Setting  this  equation  equal  to  zero  defines  the  temperature  at  which 
the  maximum  rate  occurs. 

M_  =  AnCL-xjS-1  e^a^m  (11) 

1^  is  the  sample  temperature  at  which  the  peak  in  the  DTA  thermogram 
occurs.  The  quantity  of  material  left  unreacted,  (l-x)m,  is  not  de- 
termined by  the  DTA  method. 

To  obtain  the  extent  of  reaction  as  a  function  of  temperature, 
Equation  9  can  be  integrated.  The  resulting  exponential  integral  does 
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not  yield  a  simple  expression.  Previous  work  has  resulted  in  a  satis- 
factory approximation  by  successive  integration  by  parts.  Higher  order 
terms  were  neglected,  and  the  simplified  expression  becomes; 

n(l-x)S_1  =  1  +  (n-1)  ^RTm  (12) 

^a 

The  heating  rate  does  not  appear  except  as  a  result  of  Tm,  which  varies 

with  heating  rate.  The  product  n(l-x)$~-'-  is  independent  of  heating  rate 

and  nearly  equal  to  unity.  Equation  11  now  becomes; 

0  =  ART^  e-Ea/RTm  (13) 

Ea 

Taking  the  natural  logarithm  of  both  sides  and  differentiating; 

d  ln(4>/Tm)   d  In  (AR/Ea  e-Ea/RTm) 
d(l/Tm)        d(l/Tm) 


W 


ar;      veo/  \  r 


-a 


d  InU/ffl   "Ea 
d(1000/Tm)  "  R 


By  experimentally  obtaining  DTA  thermograms  at  different  heating  rates, 
the  activation  energy  for  the  decomposition  reaction  can  be  determined 
by  Equation  14,  regardless  of  reaction  order. 
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IV.  EXPERIMENTAL  PROCEDURE 

A.  APPARATUS 

The  differential  thermal  analysis  reaction  was  obtained  by  use  of 
the  DuPont  900  Differential  Thermal  Analyzer  shown  in  Figure  1.  A  re- 
mote cell  unit,  shown  in  Figure  2,  was  used  in  this  work  due  to  the 
nature  of  the  samples  being  studied. 

DTA  uses  a  differential  thermocouple  arrangement  consisting  of  two 
thermocouples  wired  in  opposition  as  shown  in  Figure  3« 

+ 


SAMPLE 


REFERENCE       (  «Qj[ 


AT 


+ 

FIGURE  3  DIFFERENTIAL  THERMOCOUPLE 
Thermocouple  A  is  placed  in  the  sample  to  be  analyzed.  Thermocouple  B 
is  placed  in  an  inert  reference  material,  when  the  temperature  of  the 
sample  equals  the  temperature  of  the  reference,  the  thermocouples  pro- 
duce identical  voltage,  and  the  net  voltage  output,  AT,  is  zero,  when 
a  physical  or  chemical  reaction  occurs  in  the  sample,  a  differential 
signal  is  received  and  recorded.  Figure  4  shows  a  schematic  of  the  DTA 
system. 

A  silver  heating  block  is  mounted  in  the  remote  cell  assembly. 
This  block  contains  the  sample,  reference,  and  controlling  chromel- 
alumel  thermocouples  as  shown  in  Figure  5-  The  block  is  designed  to 
provide  for  equal  heating  rates  for  both  the  reference  and  sample  ma- 
terials. The  temperature  of  the  reference  is  less  than  the  temperature 
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FIGURE  2:  Remote  Cell  Assembly 
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FIGURE  4 :  DTA  System  Schematic 
(taken  from  ref .  10,  pg.  3-2) 
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FIGURE  5:  Heating  Block  and  Thermocouple  Arrangement 
(taken  from  ref.  10,  pg.  3-4) 
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of  the  block  by  a  constant  amount  that  is  directly  proportional  to  the 
heating  rate,  weight  and  heat  capacity  of  the  reference  material,  and 
inversely  proportional  to  the  thermal  conductivity  of  the  reference  ma- 
terial, plus  a  complex  exponential  function  which  can  be  neglected  for 
simplification. 

Vtb-*r(Wc/k)r 

Assuming  a  thermally  inert  sample, 

Ts=TB-<j,s(Wc/<)s 

T  =  temperature  <  =  thermal  conductivity 

4>  =  heating  rate  S  =  sample 

W  =  weight  R  =  reference 

c  =  heat  capacity  B  =  block 

The  differential  temperature  is: 

AT=TS-TR*<J>R  (Wc/k  )R-<f>3  (Wc/k  )  s 
The  heating  block  in  the  DuPont  900  is  designed  such  that  <j>q=4>r: 

AT* (wc/k )R- (Wc/k )s 
If  the  weight,  heat  capacity  and  thermal  conductivity  of  the  sample  and 
reference  material  are  the  same,  AT  is  zero.  The  AT  signal  as  a  func- 
tion of  sample  temperature  is  sent  to  the  X-Y  recorder.  Any  change  in 
the  properties  of  the  sample  results  in  a  positive  or  negative  value 
for  AT  as  shown  in  Figure  6. 


POSITIVE  SIGNAL:   (Wc/k )r> (Wc/k )s 

c,K^f(T) 
ZERO  SIGNAL:   (Wc/k )r= (Wc/k )g 

NEGATIVE  SIGNAL:   (Wc/k )r< (Wc/k )s 

C,K^f(T) 


SAMPLE  TEMP.  + 

FIGURE  6  DTA  THERMOGRAM 
21 


A  thermally  active  material  will  experience  a  change  in  heat  capac- 
ity and  thermal  conductivity  over  a  temperature  range.  The  AT  trace 
will  be  a  sloping  straight  line  if  the  properties  change  linearly  with 
temperature,  and  curved  if  the  change  with  temperature  is  quadratic  or 
of  higher  order.  Figure  7  exhibits  this  relationship,  where  the  curved 
line  would  have  a  positive  slope  if  (Wc/k)r>(Wc/k)s. 
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(Wc/k)r<(Wc/k)s 
_c,k=  linear  f(T) 
'  (Wc/k)r<(Wc/k)s 
c,k=  quadratic  f(T) 


SAMPLE  TEMP  -> 

FIGURE  7  DTA  THERMOGRAM 
The  DTA  thermogram  of  a  thermally  active  sample  will  exhibit  endo- 
therms  and  exotherms  where  the  temperatures  of  interest  are  defined  in 
Figure  8. 


It 


ONSET 


l  EXTRAPOLATED  ONSET 


SAMPLE  TEMPERATURE  + 
FIGURE  8  TRANSITION  TEMPERATURES 
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TRANSITION  TEMPERATURE  DEFINITIONS  USED  IN  FIGURE  8. 

ONSET:  The  temperature  at  which  the  thermogram  starts  to  depart  from 
the  base  line. 

EXTRAPOLATED  ONSET:  The  temperature  corresponding  to  the  intersection 

of  extrapolations  of  the  base  line  and  the  longest  straight-line 
section  of  the  low-temperature  side  of  the  peak. 

PEAK:   The  temperature  of  reversal. 

RECOVERY :  The  temperature  at  which  the  thermogram  returns  to  either 
the  same  or  a  different  base  line. 

EQUTT.TRRIUM:  The  region  where  the  state  or  form  of  the  sample  before 
transition  is  in  equilibrium  with  its  state  or  form  after 
transition. 

For  the  purposes  of  this  work,  a  fourth  well  was  drilled  in  the 
silver  heating  blocks  to  cross  check  the  temperature  scale  on  the  DTA 
900  X-Y  plotter.  A  single  chromel-alumel  thermocouple  was  set  in  inert 
glass  beads  and  connected  to  a  Varian  recorder.  This  gave  a  millivolt 
versus  time  chart  to  check  the  heating  block  temperatures  and  heating 
rate.  The  Varian  thermocouple  was  calibrated  with  a  millivolt  poten- 
tiometer and  temperatures  were  verified  through  the  use  of  a  hot  oil 
bath.  The  temperature  scale  of  the  DTA  thermogram  is  calibrated  to 
the  chromel-alumel  millivolt  output  to  give  direct  temperature  readings 
in  degrees  centigrade. 

B.  SAMPLE  PREPARATION 

The  samples  of  RDX  that  were  used  in  the  experimental  runs  were 
prepared  by  drying  a  small  quantity  of  pure  RDX  that  was  stored  in 
water.  The  sample  appeared  to  be  pure  RDX  as  its  melting  temperature 
of  204°C  compared  favorably  with  listed  values. 

Two  separate  heating  blocks  were  used  in  the  experimental  work.  In 
-the  variable  heating  rate  runs,  a  silver  block  with  wells  to  accept  two 
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millimeter  diameter  capillary  tubes  was  used.  The  RDX  sample  was 
tapped  down  into  one  of  the  four  capillary  tubes  to  a  depth  of  four 
millimeters,  which  was  calculated  to  be  19.3  milligrams  of  RDX.  The 
three  remaining  capillary  tubes  contained  equal  levels  of  inert  finely 
ground  glass  beads  .05  millimeters  in  diameter.  These  tubes  contained 
the  reference,  controlling,  and  Varian  recorder  chromel-alumel  thermo- 
couples. All  four  capillary  tubes  were  then  cut  to  1-1/8  inch.  Care 
was  taken  when  placing  the  thermocouples  to  insure  that  they  were 
pushed  to  the  full  depth  of  the  material.  The  capillary  tubes  were 
then  placed  in  their  respective  wells,  and  the  thermocouple  leads  were 
folded  under  the  rim  of  the  remote  cell  unit  to  insure  that  they  re- 
mained completely  seated  throughout  the  run.  This  was  especially  im- 
portant in  the  case  of  the  RDX  sample  since  the  reactant  gases  given 
off  during  the  exothermic  decomposition  tended  to  push  the  thermo- 
couple out  of  the  tube.  The  sample  depth  of  four  millimeters  was 
critical  for  satisfactory  results.  A  quantity  larger  than  this  would 
push  the  thermocouple  out  of  the  tube,  giving  erroneous  readings.  A 
smaller  quantity  would  not  provide  complete  coverage  of  the  thermo- 
couple probe.  The  fineness  or  texture  of  the  sample  in  this  work  was 
not  considered  since  the  reactions  occurred  above  the  melting  tempera- 
ture, or  in  all  cases  for  the  variable  heating  rate  runs,  in  the  liquid 
phase. 

The  second  silver  heating  block  was  identical  to  the  first  with  the 
exception  of  the  sample  and  reference  wells  which  were  designed  to  hold 
four  millimeter  diameter  test  tubes.  The  RDX  sample  was  measured  in 
the  same  manner  as  before  to  a  depth  of  four  millimeters,  which  in 
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this  case  gave  a  sample  of  43.32  milligrams  of  RDX.  The  second  four 
millimeter  test  tube  contained  an  equal  volume  of  inert  glass  beads  and 
served  as  the  reference.  Both  four  millimeter  tubes  were  covered  with 
small  metal  caps  through  which  the  thermocouples  were  inserted.  This 
reduced  losses  to  the  atmosphere  during  reaction.  The  controlling  and 
Varian  thermocouples  were  inserted  in  the  remaining  two  small  capillary 
wells  in  the  same  manner  as  in  the  variable  heating  rate  block.  This 
second  heating  block  was  used  in  the  isothermal  portion  of  the 
experimental  work. 

C.  VARIABLE  HEATING  RATE  EXPERIMENT 

The  DTA  900  control  panel  contains  dials  to  set  starting  tempera- 
tures and  heating  rates.  It  was  determined  that  there  was  no  apparent 
difference  in  the  thermogram  presentation  for  a  sample  that  was  heated 
from  room  temperature  at  an  established  rate  through  decomposition  or 
for  one  which  was  heated  at  an  accelerated  rate  (50°C/min.)  to  200°C, 
then  slowed  to  the  established  rate  through  the  decomposition  reaction. 
In  order  to  expedite  the  runs,  this  latter  procedure  was  adopted. 

The  full  range  of  the  DTA  900  heating  rate  capability  was  utilized. 
In  order  to  establish  consistency  in  the  procedure ,  an  outline  was  set 
down  and  followed  in  each  run. 

EXPERIMENTAL  PROCEDURE  FOR  VARIABLE  HEATING  RATE  EXPERIMENT 

(1)  Turn  DTA  900  power  knob  to  STANDBY,  allow  five  minute  warm  up. 

(2)  Place  inert  reference  and  sample  material  in  the  capillary  tubes, 
tap  and  measure  to  a  depth  of  four  millimeters.  Cut  to  1-1/8  inch. 

(3)  Insert  corresponding  thermocouples,  insuring  that  the  thermocouple 
probe  is  completely  immersed  in  the  material.  A  new  thermocouple 
is  used  for  the  RDX  sample  in  each  run. 
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(4)  Insert  capillary  tubes  into  the  correct  wells  of  the  heating  block. 

(5)  Bend  the  leads  of  the  thermocouples  under  the  remote  cell  rim. 

(6)  Insure  that  the  zero  temperature  ice  bath  reference  thermocouples 
are  fully  inserted  into  their  glass  tubes,  which  in  turn  are  im- 
mersed in  the  ice  and  distilled  water  solution  at  0°C. 

(7)  Place  a  thermogram  grid  on  the  X-Y  plotter. 

(8)  Set  the  pen  to  room  temperature  on  the  T-scale. 

(9)  Set  the  DTA  900  upper  panel  knobs: 

(a)  T-zero  shift. 

(b)  T-scale,  °C/in. 

(c)  AT-zero  shift. 

(d)  AT-scale,  °C/in. 

(e)  Baseline  slope. 

(10)  Set  the  DTA  900  lower  panel  dials: 

(a)  Temperature  heating  rate,  range  0-30°C/min. 

(b)  Set  starting  temperature  dial  to  minimum  setting  to  get  the 
nn\nimum  voltage  reading  on  the  heater  voltage  indicator. 

(11)  Set  the  program  mode  to  HEAT. 

(12)  Push  RESET  to  erase  previous  program  memory. 

(13)  Place  recorder  pen  to  the  DOWN  position. 

(14)  Place  power  knob  to  RECORD. 

(15)  Start  timer  and  Varian  recorder. 

Figure  9  shows  a  sample  thermogram  for  a  heating  rate  of  20°C/min. 
beginning  at  room  temperature.  The  endotherm  at  204°C  shows  transfor- 
mation to  the  liquid  phase,  followed  immediately  by  the  exothermic  de- 
composition with  the  peak  in  the  exotherm  at  251°C.  Figure  10  shows  a 
thermogram  for  ammonium  nitrate.  The  endotherm  at  40°C  shows  a  transi- 
tion from  rhombic  I  to  rhombic  II.  At  91°C,  a  crystalline  change  to 
tetragonal  occurs,  then  to  cubic  at  130°C.  Melting  occurs  at  174°C 
followed  by  decomposition  in  the  liquid  phase.  This  run  was  made  to 
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check  the  calibration  of  the  machine,  to  give  an  indication  of  the 
sensitivity  of  the  analyzer,  and  for  use  in  a  method  of  estimating  the 
heat  of  reaction  of  RDX,  to  be  described  in  the  recommendations.  The 
thermogram  was  compared  with  the  ammonium  nitrate  thermogram  presented 
on  page  5-7  of  the  DTA  Instruction  Manual  [Ref.  10].  The  comparison 
showed  a  disparity  in  temperature  for  the  peak  exotherm,  where  that  of 
Figure  10  showed  a  temperature  approximately  10°C  higher  than  that  of 
the  reference  thermogram.  This  was  corroborated  by  the  Varian  recorder, 
and  subsequent  temperature  readings  were  taken  from  the  Varian  recorder. 

D.  ISOTHERMAL  EXPERIMENT 

A  time  base  adaptor  for  the  DTA  900  provides  for  maintaining  the 
heating  block  at  set  temperatures  over  long  periods  of  time.  The  X-Y 
plotter  now  provides  for  time  in  minutes  on  the  abscissa  and  as  before, 
AT  on  the  ordinant.  Experimental  runs  were  made  using  the  four  mil- 
limeter test  tube  macro  heating  block.  A  variation  of  a  method  used 
in  determining  the  effects  of  thermal  aging  of  propellants  was  followed 
[Kuletz  and  Pakulak  1962].  This  technique  is  known  as  rrticroiinmersion 
autocatalytic  analysis  (MAA).  In  this  work,  the  sample  of  RDX  was 
dropped  into  the  well  of  the  macro  heating  block  after  the  block  had 
been  established  at  a  desired  temperature.  The  procedure  followed  is 
listed  below. 

EXPERIMENTAL  PROCEDURE  FOR  ISOTHERMAL  EXPERIMENT 

(1)  Turn  DTA  900  power  knob  to  STANDBY,  allow  one  to  two  hours  warm  up. 

(2)  Place  inert  reference  material  in  both  4mm  test  tubes,  tap  and 
measure  to  a  depth  of  four  millimeters. 

(3)  Place  metal  thermocouple  holder  caps  over  test  tubes  and  insert 
thermocouples,  insuring  complete  immersion. 
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(4)  Insert  reference  and  sample  4mm  test  tubes  in  the  large  wells  of 
the  macro  heating  block.  Two  small  capillary  tubes  are  used  for 
the  controlling  and  Varian  recorder  thermocouples. 

(5)  Bend  the  thermocouple  leads  under  the  remote  cell  rim. 

(6)  Check  the  zero  temperature  ice  bath  and  thermocouples. 

(7)  Place  a  thermogram  grid  on  the  X-Y  plotter. 

(8)  Set  the  time  base  adaptor  to  TIME. 

(9)  Set  the  DTA  900  upper  panel  knobs: 

(a)  T-zero  shift. 

(b)  T-scale,  Min./in. 

(c)  AT-zero  shift. 

(d)  AT-scale. 

(e)  Baseline  slope. 

(10)  Set  the  DTA  900  lower  panel  dials: 

(a)  Set  temperature  heating  rate  at  zero. 

(b)  Set  the  starting  temperature  dial  according  to  the  scale  of 
Figure  3G  of  section  3,  DTA  Instruction  Manual  [Ref.  10]. 

(11)  Set  the  program  mode  to  ISOTHERMAL. 

(12)  Allow  the  temperature  to  stabilize  at  the  desired  setting  and 
record  the  Varian  reading. 

(13)  Place  a  sample  test  tube  of  RDX  in  the  macro  heating  block  and 
reset  the  corresponding  thermocouple.   (Wear  protective  shield) 

(14)  Push  RESET  to  start  zero  tijme  at  the  left  end  of  the  scale. 

(15)  Place  recorder  pen  to  the  DOWN  position. 

(16)  Record  time  and  start  timer. 

(17)  Monitor  Varian  recorder. 

The  resulting  thermogram  plots  AT  versus  time. 
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V.  RESULTS  AMD  DISCUSSION 

The  series  of  thermograms  obtained  from  the  variable  heating  rate 
method  using  RDX  samples  are  reproduced  In  Figures  11-17.  During  each 
of  these  runs,  temperatures  were  recorded  at  5°C  intervals  and  at  sig- 
nificant points  on  the  thermogram  (melting,  peak  exotherm)  with  the 
corresponding  time  in  seconds.  This  data  is  plotted  in  Figures  18-36. 
An  n^*1  order  regression  analysis  program  of  the  WANG  700  was  used  to 
calculate  the  least  squares  fit  to  the  data  points,  the  slope  of  which 
is  the  actual  heating  rate  for  that  run.  The  resulting  equation  gives 
the  temperature  at  any  time  during  each  run  for  the  programed  heating 
rate.  The  calculated  heating  rates,  along  with  the  sample  temperature 
of  the  peak  in  the  exotherm  Tm,  are  the  required  data  for  determining 
the  kinetic  parameters  of  the  decomposition  by  Kissinger's  method  [Ref. 
2].  These  data  are  summarized  in  Table  I. 

The  kinetic  data  for  the  decomposition  reaction  of  RDX  were  deter- 
mined from  the  plot  of  the  variable  heating  rate  experimental  data  of 
Table  I.  The  natural  logarithm  of  the  rate  factor  <}>/Tm  was  plotted 
versus  the  reciprocal  of  the  Kelvin  temperature  of  the  heating  block, 
as  shown  in  Figure  37.  The  WANG  700  was  used  to  calculate  a  least 
squares  fit  to  the  experimental  data  points.  The  slope  of  this  line 
determines  the  activation  energy  which  was  calculated  to  be  49.^713 
kcal/mole.  The  Arrhenius  frequency  factor  was  determined  from  Equation 
13,  averaged  over  the  temperature  range  of  the  experiment.  The  average 
value  was  6.9199  x  1020,  which  varied  from  the  fastest  heating  rate 
value  by  1.1501  x  1020  and  from  the  lowest  by  0.1415  x  1020.  These 
values  compare  favorably  with  those  of  previous  works  [Rauch  and 
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Wainright  1969,  Robertson  1948].  The  kinetic  parameters  are  summarized 
in  Table  II.  There  is  a  1.6%   error  in  the  value  for  the  activation  en- 
ergy when  compared  with  that  of  Rauch  and  Wainright's  work  for  Picatinny 
Arsenal. 

Appendix  A  is  an  IBM  360  computer  readout  which  tabulates  values  for 
the  heating  rate  with  incremental  increase  in  the  peak  temperature  Tm.  The 
experimentally  determined  values  of  activation  energy  and  Arrhenius  fre- 
quency factor  were  used  in  Equation  13  to  generate  values  of  the  heating 
rate  <}>.  Appendix  B  tabulates  the  same  information  for  comparison  using 
the  kinetic  data  determined  by  Rauch  and  Wainright.  This  data  is  presented 
in  Table  III  along  with  values  from  the  experimental  runs  for  comparison. 
At  higher  heating  rates,  the  experimental  and  calculated  data  differ  to  a 
larger  extent.  This  is  probably  due  to  experimental  time  lag  at  the  higher 
heating  rates  which  caused  lower  temperature  readings  at  the  higher  rates. 
Although  the  values  of  activation  energies  compare  favorably,  Table  III 
shows  that  the  experimentally  determined  temperatures  of  the  peak  in  the 
exotherms  at  the  various  heating  rates  would  have  to  be  approximately  30°C 
higher  to  give  heating  rates  comparable  to  the  Picatinny  data. 

The  rate  constant  k  over  the  experimental  temperature  range  was  calcu- 
lated using  the  Arrhenius  equation.  Appendix  C  tabulates  values  of  k  using 
the  experimentally  determined  values  of  Ea  and  A.  Appendix  D  presents  the 
same  data  using  the  kinetic  data  determined  in  the  study  made  for  Picatinny 
Arsenal.  These  values  are  compared  in  Table  IV.  The  smaller  values  for 
the  rate  constant  resulting  from  the  Picatinny  data  correspond  with  the 
lower  heating  rate  values  shown  in  Table  III. 
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Sane  variation  existed  in  the  peak  temperature  of  the  exotherms  as 
evidenced  in  Table  I.  This  resulted  in  the  scatter  noted  in  Figure  37. 
An  attempt  to  reproduce  the  peak  temperature  at  a  given  heating  rate 
was  made  over  three  consecutive  runs  using  the  standard  experimental 
technique.  Results  are  shown  in  Figure  38,  where  the  three  runs  were 
carried  out  at  a  heating  rate  of  6°C/min.  The  maximum  variation  in  the 
peaks  was  4°C.  This  compares  with  the  variation  in  temperatures  for 
the  experimental  and  calculated  values  of  Tm  shown  in  Table  III.  The 
maximum  deviation  approaches  .8%   error,  which  is  well  within  the  accur- 
of  the  experiment.  There  is  a  definite  increase  in  Tm  with  heating 
rate  as  proposed  by  the  DTA  kinetic  theory  [Kissinger  1957].  The  ef- 
fect of  sample  size  appeared  to  be  only  an  increase  in  the  quantity  of 
heat  liberated,  as  shown  in  Figure  39.  The  peak  temperature  did  not 
vary  to  any  greater  degree  than  in  the  normal  manner  with  equivalent 
volumes  of  material. 

The  isothermal  data  gave  somewhat  less  conclusive  results.  The 
thermograms,  Figures  40-44,  show  a  difinite  exothermic  decomposition 
reaction  which  occurs  at  a  constant  block  temperature  after  an  induc- 
tion time  which  decreases  with  increased  temperature.  A  temperature 
range  of  192-210°C  was  investigated.  The  time  elapsed  to  initiation 
of  the  exothermic  reaction,  or  induction  time,  and  the  time  elapsed  to 
the  peak  in  the  exotherm,  the  time-to-deflagration,  were  observed  and 
recorded.  Table  V  summarizes  the  isothermal  data.  It  was  noted  that 
although  these  reactions  occur  at  or  below  the  normal  melting  tempera- 
ture they  occurred  in  the  liquid  phase.  A  physical  observation  during 
one  of  the  runs  confirmed  a  liquid  phase  decomposition  with  decomposi- 
tion gases  bubbling  off.  The  sample  at  192°C  showed  no  evidence  of 
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reaction  after  six  and  one-half  hours.  Near  the  melting  temperature 
(204°C)  and  above,  the  sample  changed  to  the  liquid  phase  immediately 
upon  reaching  block  temperature  with  the  exothermic  decomposition  en- 
suing. The  micro immersion  autocatalytic  analysis  (MAA)  method  was  used 
to  determine  kinetic  data  from  the  isothermal  phase  of  the  experiment. 
This  method  is  based  on  the  equation; 


T2  =  ^a   Ea/RT  (15) 


cR  e 

where  Q  is  the  heat  of  reaction.  A  plot  of  the  natural  logarithm  of 
the  rate  factor  (x/T2)  versus  the  reciprocal  of  the  constant  block  tem- 
perature gives  a  slope  equal  to  Ea/R.  In  addition,  an  Arrhenius  plot 
using  the  reciprocal  of  the  time  data  versus  1000/T  was  made  [Kuletz 
and  Pakulak  1962].  Both  induction  time  and  time-to-deflagration  data 
were  used.  Table  VT  lists  the  induction  time  data.  Two  plots  were 
made  using  these  data.  The  natural  logarithm  of  the  reciprocal  of  the 
induction  time  versus  the  reciprocal  of  the  isothermal  block  tempera- 
ture in  degrees  Kelvin  is  presented  in  Figure  45.  The  MAA  procedure, 
plotting  the  natural  logarithm  of  the  rate  factor  (t0/T2)  versus  1000/ 
T°K,  is  shown  in  Figure  46.  This  was  done  again  using  the  time-to- 
deflagration.  The  data  for  time-to-deflagration  are  listed  in  Table 
VTI.  The  plots  of  these  data  are  shown  in  Figures  47  and  48  respec- 
tively. The  results  are  shown  in  Table  VTII.  In  all  cases  the  WANG 
700  was  used  to  determine  a  least  squares  fit  to  the  data  points  and 
the  resulting  slopes.  It  can  be  seen  that  both  the  reciprocal  time 
and  MAA  rate  factor  methods  produce  nearly  identical  results.  The 
values  obtained  for  Ea  are  far  too  high  for  activation  energies. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 

Experimental  work  done  on  HMX  [Suryamarayana  and  Gray bush  19 67] 
suggested  that  materials  which  exhibit  sigmoid  pressure-time  curves 
generally  involve  three  stages  in  the  decomposition  reaction.  The 
first  stage  is  the  induction  period.  This  is  followed  by  an  acceler- 
ation period  where  the  rate  reaches  a  maximum..  The  material  is  con- 
sumed in  the  decay  period.  Even  with  this  mass  spectrometric  study, 
the  conclusions  were  limited  due  to  the  several  decomposition  products 
formed.  It  was  observed  that  each  product  reached  a  constant  rate  of 
formation  after  an  acceleratory  stage  which  appeared  due  to  autocat- 
alysis.  It  was  found  that  HMX  exhibited  different  modes  of  decomposi- 
tion. The  kinetic  data  changed  in  different  discrete  temperature 
ranges.  These  findings  were  supported  in  a  DTA  analysis  of  HMX  [Hondee 
1971]j  where  three  separate  activation  energies  were  found  during  three 
discrete  stages  of  the  decomposition  process. 

Causes  for  acceleration  of  the  decomposition  of  HMX  to  a  constant 
rate  can  be  attributed  to  [Suryamarayana  and  Graybush  1967]: 

(1)  Progressive  melting  as  a  result  of  lowered  melting  point  by 
the  decomposition  products. 

(2)  Self-heating. 

(3)  Autocatalysis  by  products. 

(4)  Acceleration  due  to  structural  factors  such  as  the  increase 
in  the  number  of  nuclei. 

During  this  experimental  study  of  the  decomposition  of  RDX  there 

was  no  evidence  of  different  activation  energies  since  the  heating  rate 

in  all  runs  remained  constant  over  the  entire  decomposition.  The  peak 

exotherm,  although  somewhat  variable,  did  not  exhibit  more  than  one 

sharp  peak.  The  thermograms  for  the  variable  heating  rate  experiments 
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(Figures  11-17)  show  that  the  exothermic  decomposition  reached  a  peak 
which  was  in  most  cases  followed  by  a  second  peak.  This  second  peak 
is  probably  a  result  of  the  capillary  action  of  the  liquid  sample 
causing  the  thermocouple  to  rise  in  the  tube,  then  drop  back  down  just 
before  the  reaction  completes  causing  another  jump  in  the  thermogram. 

The  isothermal  experiment  thermograms  (Figures  40-44)  show  a  def- 
inite decomposition  reaction  after  an  induction  time  at  lower  tempera- 
tures. This  is  in  agreement  with  previous -  studies  which  have  found 
evidence  of  self-heating  and  autocatalysis  in  the  decomposition  of  RDX 
and  other  explosives.  These  data  can  be  applied  to  estimates  of  the 
sensitivity  of  explosives  and  propellants. 

Several  reasons  may  be  proposed  for  the  discontinuities  in  the 
variable  heating  rate  experiment.  Although  the  RDX  sample  was  assumed 
homogeneous  and  pure,  it  may  well  have  been  nonuniform  in  nature  and 
composition.  Possible  water  content  in  the  samples  was  neglected. 
DTA  results  vary  between  laboratories  due  to  the  difference  in  experi- 
mental techniques  and  equipment  as  well  as  in  the  interpretation  of 
the  results.  Standardization  of  method  would  be  beneficial. 

Explosive  decomposition  rates  vary  rapidly  with  temperature  and 
for  ideal  results  the  temperature  of  the  sample  should  at  all  times 
be  uniform.  This  is  not  the  true  situation  due  to  thermal  conduction 
in  the  sample  and  self-cooling  or  self-heating  effects  arising  from 
the  reaction.  Robertson  [Ref.  9]  showed  that  errors  in  the  values  of 
the  kinetic  parameters  occur  due  partly  to  autocatalysis  and  self- 
heating.  Autocatalysis  increases  the  activation  energy  value  by  ap- 
proximately 2000  cal/mole  as  a  maximum  error  and  causes  a  corresponding 
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variation  in  the  Arrhenius  frequency  factor  by  as  much  as  a  factor  of 
10.  Errors  due  to  self -heating  are  also  incurred  and  tend  to  increase 
Eg^  and  A  as  well.  This  effect  increases  with  increased  temperature. 

Experimental  errors  are  also  induced  by  the  apparatus.  The  accu- 
racy of  the  DTA  method  depends  on  the  precision  with  which  a  uniform 
heating  rate  can  be  maintained.  Although  the  DTA  900  provides  for  a 
constant  rate  to  the  heating  block,  the  exterior  of  the  sample  is  hot- 
ter than  the  interior  and  reacts  sooner  and  faster.  This  causes  addi- 
tional autocatalytic  and  self-heating  effects  which  are  accented  by  an 
increased  heating  rate.  Another  source  of  error  in  this  work  was  the 
uncontrolled  atmosphere.  This  provides  for  possible  variations  in  the 
data  from  day  to  day.  Bubbles  of  the  decomposition  gases  which  formed 
on  the  thermocouples  caused  the  high  intensity  chatter  on  the  thermo- 
grams. Capillary  action  of  the  melt  during  decomposition  caused  it 
to  climb  the  walls  of  the  tube  surrounding  the  thermocouple  sleeve. 
This  drew  the  material  away  from  the  thermocouple  probe.  Care  had  to 
be  taken  to  keep  the  level  of  the  liquid  below  the  sleeve  in  order  to 
miriimize  this  effect. 

Notwithstanding  these  adverse  and  random  effects,  the  dominant  fac- 
tor in  controlling  the  shape  and  position  of  the  endotherms  and  exo- 
therms  is  the  nature  of  the  reaction  itself.  The  favorable  results 
obtained  lend  credibility  to  the  DTA  method  and  suggest  that  further 
development  in  the  technique  will  prove  useful  in  providing  a  fast  and 
reasonably  accurate  method  of  characterizing  the  thermal  and  kinetic 
properties  of  new  explosives  and  propellants. 

Areas  for  further  study  are  wide  and  varied.  Considerably  more 
information  as  to  the  mechanisms  of  the  reaction  can  be  obtained  by 
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combining  the  DTA  method  with  one  or  more  of  several  more  precise  an- 
alytical techniques  such  as  thermogravimetry ,  infrared  spectrometry, 
mass  spectrometry,  X-ray  measurements,  and  visual  observations. 

More  information  may  be  determined  from  the  DTA  thermograms  them- 
selves. It  has  been  shown  that  the  area  under  the  peak  in  the  thermo- 
gram is  proportional  to  the  heat  of  reaction  [Rivette  and  Besser  1961]. 
This  provides  a  means  to  quickly  estimate  the  heats  of  reaction  of  the 
sample.  The  shape  of  the  exotherm  provides  a  means  to  estimate  the 
reaction  order  [Kissinger  1957].  Rough  estimations  using  this  method 
showed  the  reaction  order  to  vary  from  0.2  to  0.8  over  the  range  of 
heating  rates  considered.  Further  analysis  of  the  isothermal  thermo- 
grams appears  warranted  to  discover  the  significance  of  the  two  slopes 
on  the  low  temperature  side  of  the  exotherm  and  the  shape  and  extent 
of  the  decay  pattern  on  the  high  temperature  side,  as  well  as  the  sig- 
nificance of  the  high  values  obtained  for  Ea  from  the  Arrhenius  and  • 
MAA  plots  of  the  induction  time  and  time-to-deflagration  data. 
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TABLE  I 
EXPERIMENTAL  DATA,  VARIABLE  HEATING  RATE  THERMAL  ANALYSIS 

RUN  #  <j),0C/MIN 

1  2.431 

2  2.880 

3  3-333 

4  3.8571 

5  4.3636 

6  4.6957 

7  5.1940 

8  5.6552 

9  6.2069 

10  6.8182 

11  7.9149 

12  8.8696 

13  9.7143 

14  10.5556 

15  11.6471 

16  16.5217 

17  21.0000 

18  26.0769 

19  30.0000 


Tm  °K 

2 

lnUAn) 

(1/Tm)xl03 

504.00 

-II.56OI 

1.9802 

501.00 

-11.3754 

1.9960 

509.00 

-11.2610 

1.9646 

513.00 

-11.1306 

1.9493 

505.25 

-10.9768 

1.9792 

509.00 

-10.9183 

1.9646 

510.50 

-10.8233 

1.9589 

514.50 

-10.7538 

1.9436 

515.25 

-10.6636 

1.9408 

513.00 

-10.5610 

1.9493 

517.25 

-10.4283 

1.9333 

510.75 

-10.2891 

1.9579 

511.50 

-10.2011 

1.9550 

512.75 

-10.1229 

1.9503 

519.00 

-10.0488 

1.9268 

517.00 

-9-6914 

1.9342 

524.50 

-9.4804 

1.9066 

523.25 

-9.2591 

1.9111 

528.00 

-9.1370 

1.8939 
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TABLE  II 

KINETIC  PARAMETERS 

Ea(kcal/mole) 

EXPERIMENTAL 

49-5 

RAUCH  &  WAINRIGHT 

48.7 

ROBERTSON 

47.5 

A  (sec-1) 

6.9X1020  1020.8 
1.7X1019  1019.2 
3-16X1018   1018.5 


TABLE  III 

HEATING  RATE  DATA  COMPARISON 


PEAK 
TEMPERATURE 

HEATING  RATE 

* 

Tm 

EXPERIMENTAL  CALCULATE 
2.4310      2.4817 

D  PICATINNY 

504.00 

0.1338 

509.00 

4.6957      4.1122 

0.2200 

513.00 

6.8182      6.1167 

0.3253 

524.50 

21.0000     18.5317 

0.9694 

528.00 

30.0000     25.7247 
TABLE  IV 

1.3391 

RATE  CONSTANT  DATA  COMPARISON 

PEAK 
TEMPERATURE     RATE  CONSTANT  k 

(sec  ) 

Tm 

EXPERIMENTAL    I 
0.2433 

ICATINNY 

504.00 

0.0129 

509.00 

0.3952 

0.0208 

513.00 

0.5787 

0.0303 

524.50 

1.6772 

0.0864 

528.00 

2.2974 

0.1177 
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TABLE  V 

EXPERIMENTAL  DATA,  ISOTHERMAL  ANALYSIS 

INDUCTION  TIME     DEFLAGRATION 
RUN  #      T  (°C)   T  (°K)   t0  (min)  t0  (sec)   t  (min)  t  (sec)   ATmax 


1 

192.00 

465.00 

EXCESS  OF  6-1/2  HRS. 

2 

194.25 

467.25 

192.0 

11,520 

238.5 

14,310 

.39 

3 

196.75 

469.75 

130.0 

7,848 

145.6 

8,736 

.30 

4 

198.00 

471.00 

85.0 

5,100 

107.5 

6,450 

.66 

5 

198.50 

471.50 

37-5 

2,250 

58.5 

3,510 

1.30 

6 

199.50 

472.50 

15.0 

900 

34.8 

2,088 

1.50 

7 

204.50 

477.50 

10.8 

648 

19.2 

1,152 

1.60 

8 

206.75 

479.75 

0.0 

0 

0.0 

0 

3.20 

TABLE  VI 

ISOTHERMAL  ANALYSIS,  INDUCTION  TIME  DATA 

RUN  #         1000/T    t0  sec      lAo      ln(l/t0)    ln(t0/T2) 


2 

2.1402 

11,520 

.868lxl0-1+ 

-9.3518 

-2.9419 

3 

2.1288 

7,848 

. 1274x1 0"3 

-8.9680 

-3-3364 

4 

2.1231 

5,100 

. 1961x1 0"3 

-8.5370 

-3.7727 

5 

2.1209 

2,250 

.4444xlO~3 

-7.7187 

-4.5932 

6 

2.1164 

900 

.llllxlO-2 

-6.8024 

-5.5137 

7 

2.0964 

648 

. 1543x1 0~2 

-6.4739 

-5.8611 
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TABLE  VII 
ISOTHERMAL  ANALYSIS,  TIME-TO-DEFLAGRATION 


RUN  # 

1000/T 

lACsec"1) 

ln(x/T2) 

ln(l/x) 

2 

2.1402 

.6988xlO_l+ 

-6.8194 

-5.4744 

3 

2.1288 

. 1145x10" 3 

-7.3235 

-4.9809 

4 

2.1231 

. 1550x10" 3 

-7^6322 

-4.6775 

5 

2.1209 

.2849xl0"3 

-8.2428 

-4.0690 

6 

2.1164 

.4789xl0"3 

-8.7665 

-3.5496 

7 

2.0964 

. 8681x10" 3 
TART,E  VIII 

-9.3801 

-2.9549 

ISOTHERMAL  RESULTS 

INDUCTION  TIME    TIME-TO-DEFLAGRATION 
1/tp     tp/T2     1/t      t/T2 


Ea  (kcal/mole)->  159-0282  160.9012   131.0144  132.8778 
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Figure  12.  Variable  Heating  rate  runs  4  8i5. 
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Figure  13.  Variable  Heating  Rate  Runs  6-8. 
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Figure  14.  Variable  Heating  Rate  Runs  9-11. 
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Figure  15.  Variable  Heating  Rate  Runs   12-14. 
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Figure  16.  Variable  Heating  Rate   Runs  15-17. 
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APPENDIX  A 
HEATING  RATE  VS.  PEAK  TEMPERATURE,   500. 00-530. 50°K 

[EATING  RATE  =  #  e"Ea^Tm 
Ea 

A  =  6.91985  x  1020 

Ea  =  49,471.30  cal/mole 

R  =  1.987  cal/mole  °K 


TEMP.    DEG    K  RATH  TEMP.    DEG    K  RATS 


5CC.00      -  1.64513  515.50  7.81557 

500. 50      1.73252  516.00  -    8.20592 

501.00      -  1.82433  516.50  8.61496 

5C1.50      -  1.92082  517.00  9.04354 

502.  OC 2.02223  517.50 9.49271 

502.50      -  2.12876  518. CO  -    9.96309 

503.00      -  2.24069  518.50  10.45597 

503.50      2.35825  519.00  10.97204 

504. GO      -  2.48174  519.50  11.51268 

5C4.50      2.61146  520.00  -    12.07884 

505.00      2.74766  520.50  12.67184 

505.50      -  2.89069  521.00  -    13.29251 

5C6.00      -  3.04084  521.50  -    -    -  13.94250 

506.50      -  3.19853  522.00  -    14.62290 

5C7.00      -  3.36404  522.50  -    -    -  15.33508 

507.50      -  3.53773  523.00  -    16.08067 

508.00      -  3.72009  523.50  -    -    -  16.86069 

5C8.50 3.91143  524.00 17.67717 

509.00      -  A. 11222  524.50  18.53172 

509.50      -  4.32292  525.00  -    19.42580 

51C.00  -    -    -  4.54399  525.50  -    -    -  20.36110 

510.50       4.77592  526. CO  -    21.33943 

511.00      -  5.01914  526.50  -    22.36304 

511.50      -  5.27426  527.00  -  23.43390 

512.00      -  5.54182  527.50  -    -    -  24.55339 

512.50      -  5.82250  528.00  -    -    -  25.72473 

513.00      -  6.11673  528.50  -  26.94904 

513.50      -  6.42524  529.00  -    -    -  28.22940 

514.00  -    -    -  6.74866  529.50  -    -    -  29.56830 

514.50       -  7.08771  530.00  -    -    -  20.96779 

515.00 7.44310  530.50 32.43098 
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APPENDIX  B 

HEATING  RATE  VS.   PEAK  TEMPERATURE,   500. 00-530. 50°K 

TTT,_TO  „nmr.         ARTm  p-Ea/RTm 
HEATING  RATE  =     -g—  e 

A  =  1.70000  x  1019 
Ea  =  48,700.00  cal/mole 
R  =  1.987  cal/mole  °K 


TEMP.    DEG   K  RATE  TEMP.    DEG    K  RATE 


500.00  C. 08924  515.50      0.41415 

5C0.50 0.G9390  516. CO 0.43452 

5C1.00  -  0.C9880  516.50       0.45585 

501.50  0.10395  517.00      0.47818 

502.00  -  0.10935  517.50      0.50156 

502.50  -  0.11502  518.00      0.52603 

503.00  -  0.12098  518.50  -    -    -  0.55165 

5C3.5C  -  0.12723  519.00       0.57847 

504.00  -  0.13379  519.50      -    0.60654 

504.50  -  0.14067  520.00      -  0.63590 

5C5.C0  -  0.14790  520.50  -    -    -  0.66665 

505.50  -  0.15548  521.00      0.69880 

5C6.00  -  0.16343  521.50       0.73244 

506.50  -  0.17177  522.00      0.76763 

507.00  -  0.18052  522.50      0.80446 

5C7.50  -  0.18971  523.00  -    -    -  0.84296 

508.00  -  0.19933  523.50      -  0.88323 

5C8.5C  -  0.20943  524.00       -    0.92535 

5C9.CC  -  0.22002  524.50  -    -    -  0.96939 

509.50  -  0.23112  525. CO      -    1.01545 

510.00  -  0.24275  525.50      1.06359 

510.50  -  0.25495  526.00      -  1.11392 

511.00  -  0.26774  526.50  -    -    -  1.16653 

511.50  -  0.28114  527.00  -    -    -  1.22153 

512.00  -  0.29518  527.50      -  1.27901 

512.50  -  0.30990  528.00  -    -    -  1.33906 

513.00  -  0.32532  528.50  -    -    -  1.40185 

513.50  -  0.3^147  529. CC  -    -    -  1.46742 

514.00  -  0.35840  529.50  -    -    -  1.53594 

514.50  -  0.37613  530.00  -    -    -  1.60753 

515.00  -  0.39470  530.50  -    -    -  1.68233 
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APPENDIX  C 

RATE  CONSTANT  VS.   ISOTHERMAL  TEMPERATURE,   470.00-550.00  °K 

k  =  A  e-Ea/Rr 
A  =  6.91985  x  1020 

Ea  =  49,471.30  cal/mole 

R  =  1.987  cal/mole  °K 


TEMP.  05G  K 


470.00 
470.50 
471.00 
471. 5C 
472.00 
472.50 
473.00 
473.50 
474.00 
474.50 
475.0C 
475.50 
476.00 
476.50 
477.00 
477.  5C 
478.00 
478.50 
479.00 
479.50 
480.00 
48C.50 
481.00 
481.50 
482.00 
482.50 
4  8  3.  OC 
483.50 
484.00 
484.50 


C. 00682 
0.00722 
0.00764 
0.00  80  8 
0.00854 
0.00903 
0.00955 
0.01010 
0.01067 
0.01128 
0.01192 
0.01259 
0.01331 
C. 01406 
0.01485 
0.C1568 
0.01656 
0.01749 
0.01846 
0.01949 
0.02058 
C. 02172 
C. 02292 
0.02418 
0.02552 
0.02692 
0.02  840 
0.02995 
0.C3159 
0.03331 


TEMP.  OEG  K 


510.50 
511.00 
511.50 
512.00 
512.50 
513.00 
513.50 
514.00 
514.50 
515.00 
515.50 
516. CO 
516.50 
517.00 
517.50 
518.00 
518.50 
519.00 
519.50 
520. OC 
520.50 
521.00 
521.50 
522.00 
522.50 
523.00 
523.50 
524. CC 
524.50 
525.00 


0.45627 
0.47857 
0.50191 
0. 52634 
0.55192 
0.57868 
0. 60669 
0.63599 
C. 66664 
0.69871 
0.73225 
0.76733 
0.80402 
0.84239 
0. 88252 
0.92446 
0.96333 
1.01416 
1.06209 
1.11218 
1.16454 
1.21923 
1.27640 
1.33613 
1.39852 
1.46371 
1.53179 
1.60290 
1.67718 
1.75475 
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APPENDIX  C   (continued) 


TEMP.    DEG    K  K  TEMP,    DEG   K 


485.00  0.03512  525.50  1.83574 

485.50  -  0.03703  526.00  1.92029 

486.00  -  0.03904  526.50  -  2.00858 

486.50 0.04115  527.00 2.10077 

487.00 0.04337  527.50 2.19696 

487.50  -  0.04570  528.00  -    2.29741 

488.00  -  0.04816  528.50  -    -    -  2.40220 

488.50  -  0.05074  529. CO  -  2.51157 

489.00  -  0.05345  529.50  2.62573 

489.50  -  0.05630  530.00  -    2.74482 

490.00  -  0.05930  530.50  2.86910 

49G.50  -    -    -  0.06246  531.00  -    -    -  2.99872 

491.00  -  0.06577  531.50  2.13397 

491.50  -  0.06925  532.00  2.27501 

492.00  -  0.07291  532.50  -  3.42214 

492.50  -  0.07675  533. OC  3.57555 

493.00  -  0.08079           •  533.50  -    -    -  3.73556 

492.50  -  0.C8503  534.00  -    3.90243 

494.00  -  0.08949  534.50  -    -    -  4.07638 

494.50 0.09417  535. OC  -    -    -  4.25776 

495.00  0.09908  535.50  -    4.44687 

495.50  -  0. 10424  536.00  -    4.64402 

496.00  -  0.10965  536.50  -    -    -  4.84947 

496.50  0.11534  537. CC  -  5.06354 

497.00  -  0.12131  537.50  -    5.28675 

497.50  -  0.12757  538.00  5.51929 

498.00  0.13414  538.50  5.76161 

498.50  -  0.14104  539.00  6.01412 

499.00  -  0.14828  539.50  6.27721 

499.50  -  0.15537  540.00  6.55121 

500.00  -    -    -  0.16384  540.50  -    -    -  6.83666 

500.50  0.17220  541.00  7.13399 

501.00  -  0.18096  541.50  -    7.44369 

501.50  -  0.19015  542.00  -    7.76624 

502.00  -  0.19979  542.50  -  8.10215 

502.50  -  0.20990  543.00  8.45181 

503.00  0.22050  543.50  -    -    -  8.81602 

503.50  -  0.23160  544.00  -    -    -  9.19509 

5C4.CC  -  0.24325  544.50  -  9.58973 

504.50  -  0.25546  545.00  -    -    -  10.00054 

505.00  -  0.26825  545.50  10.42831 

505.50  -  C. 28165  546. CO  -    10.87339 

506.00  -  0.29570  546.50  -    -    -  11.33659 

506.50  -  0.31042  547.00  -    -    -  11.81863 

5C7.00  -  0.32584  547.50  -    -    -  12.32022 

507.50  -  0.34199  548.00  -  12.84212 

5C8.00  -  C. 35891  548.50  -  13.38532 

508.50  0.37663  549.00  -    13.95021 

5C9.00  -  0.39518  549.50  14.53783 

5C9.50  -  0.41461  550.00  -    15.14928 
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APPENDIX  D 

RATE  CONSTANT  VS.   ISOTHERMAL  TEMPERATURE,   470.00-550.00  °K 

k=Ae-Ea/RT 

A  =  1.7000  x  1019 

Ea  =  48,700.00  cal/mole 

R  =  1.987  cal/mole  °K 


TEMP.    DEG    K  K  TEMP.    DEG    K  K 


470.00      -  0.00038  510.50      0.02398 

470. 5C 0.00040  511.00  -    -    -  0.02513 

471.00      0.00043  511.50      0.02634 

471.50      -  0.0C045  512.00      0.02760 

472.00      -  0.00048  512.50  -    -    -  0.02892 

472.50      -  0.00050  513.00      0.03030 

473.00      -  0.00053  513.50       -    C. 03174 

473.50      -    0.00056  514.00      0.03325 

474.00      0.00059  514.50      0.03483 

474.50      -  0.0C063  515.00       -    0.03647 

475.00      -  0.00066  515.50  -    -    -  0.03820 

475.50      -  0.00070  516.00      0.C40C0 

476.00      -  C. 00074  516.50  -    -    -  0.04188 

476.50  -    -    -  0.00078  517. OC 0.04385 

477.00      -  0.00082  517.50       0.04590 

477.50      -  0.CO087  518.00      0.04805 

478.00 0.00092  518.50  -    -    -  C. 05029 

478.50 0.00097  519.00  -    -    -  0.05264 

479.00 0.00102  519.50 0.05508 

479.50      -  0.00108  520.00      -    0.05764 

48C.CC      -  0.00113  520.50  -    -    -  0.06031 

480.50      -  0.00120  521.00      -    0.06310 

481.00      -  0.00126  521.50       -  0.06601 

481.50  -    -    -  0.00133  522.00      -    0.06905 

482.00      -  0.00140  522.50      -  0.07222 

482.50      -  0.00148  523.00  -    -    -  0.07553 

482.00      -  0.00156  523.50      -    0.07899 

483.50      -  0.00164  524.00       -    0.08260 

484.00      -  0.C0173  524.50  -    -    -  0.08637 

484.50      0.00182  525.00      C. 09030 
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APPENDIX  D  (continued) 


TEMP.    DEG   K 


TEMP.    DEG    K 


co  c    en      -    -    -         0.094-40 

Aflf-   ftn 0.00192  112'nn 0.09868 

J|5.00      0.00202  526.C0      0.10314 

1|5.50      0.00213  l?700      0.10780 

Afi6*50      0.C0224  527 .00      0.11266 

111   00 0.00236  I?p    00 0.11772 

1?VSQ 0.00249  |28 .00 0.12301 

iU'cO 0.00262  lie    CO 0.12852 

IflP    50 0.00  276  5|9 .00 0.13427 

aaQ   00 0.00290  529 .5Q 0.14026 

Apq-°50      0.C0306  530 .00       0.14651 

ac?oo 0.00322  111    no 0.15303 

Aon    50 0.00339  531. 00 0.15982 

lq?    00 0.00356  531 .50 0.16690 

lll'Sn      0.00375  111*50      0.17427 

5qp*0  0 0.00394  III   00 0.18196 

aq?'°50 0.00415  \\V%0 0.18998 

tll'lo 0.00436  III    00 0.19833 

tli    50 0.00459  I?450 0.20703 

HI' no 0.00482  %\%ar\ 0.21609 

JlJ'Sn      0.CC507  ||l-2§      0.22553 

ac^oO 0.00533  III  00 0.23537 

ill'28      0.00561  |36 .00      c. 24562 

211*00 0.00589  III- 00 0.25629 

iIt#Y=? 0.00619  11-750 0.26741 

JS$#no      C.00651  lll'no      0.27899 

22?*92      0.0C684  |l§*Sn       0.29104 

til- CO      6.00719  llg*00      0.30359 

498    50 0.00755  III    50 0.31665 

III   00 0.00793  111    00 0.33026 

JZl'w -      0.00833  cln    FT C.  34442 

^'no 0.00875  |??#on 0.35916 

500*50 0.00919  541  .  00 0.37451 

501    00 0.CO965  caI'oO 0.39047 

50 I    50 0.01013  It?    SO 0.40710 

In?    00 0.01064  VAoo 0.42439 

18 2    50 0.01U6  111    50 0.44238 

5ol    00      0.01172  Iaa'oO       0.46110 

lnV?n 0.01230  Iaa*5P 0.48058 

!8!-00      0.01291  |ar-00       0.50084 

«C4.50 0.01355  545. 00 0.52192 

505   00 0.01421  IaaqO 0.54383 

505    50 0.01491  12a    50 0.56664 

5? ft    00 0.01564  liVnn 0.59034 

?£f  50 0.01641  Ia7    50 0.61501 

5?7*00 0.01721  HI    no 0.64065 

5C750 0.01805  Iap-5C 0.66731 

508   00 0.01893  !aq-00 0.69502 

508.50      0.01985  raq-50      0.72384 

lol'.OO 0.02081  549.50 0#75378 

5C9.5C "      O.O^ia^ 
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